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Abstract 

We numerically investigate the features of typical orbits occurring 
in the Oort cloud (r ss 50 — 150 kAU) in the low-acceleration regime 
of the MOdified Newtonian Dynamics (MOND). We take into account 
the so-called External Field Effect (EFE) because the solar system is 
embedded in the Milky Way. In the framework of MOND this does 
matter since the gravitational acceleration of Galactic origin felt by 
the solar system is of the same order of magnitude of the characteristic 
MOND acceleration scale ^4o ~ 10~ 10 m s -2 . We use three different 
forms of the MOND interpolating function fj,[x) and two different val- 
ues for the Galactic field at the Sun's location. We find that MOND 
produces highly distorted trajectories with respect to the Newtonian 
case, especially for very eccentric orbits. It turns out that the shape 
of the MOND orbits strongly depend on the initial conditions. For 
particular initial state vectors, the MOND paths in the ecliptic plane 
get shrunk extending over much smaller spatial regions than in the 
Newtonian case, and experience high frequency variations over one 
Keplerian orbital period. Ecliptic orbits with different initial condi- 
tions and nearly polar orbits are quite different getting distorted as 
well, but they occupy more extended spatial regions. This facts may 
have consequences on the composition and the dynamical history of 
the Oort cloud which are difficult to predict. 
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1 Introduction 

In many astrophysical systems like, e.g., spiral galaxies and clusters of galax- 
ies a discrepancy between the observed kinematics of their exterior parts and 
the predicted one on the basis of the Newtonian dynamics and the matter 
detected from the emitted electromagnetic radiation (visible stars and gas 
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clouds) was present since the pioneering studies by Bosnia pQ and Rubin 
and coworkers [2 J on spiral galaxies. More precisely, such an effect shows 
up in the galactic velocity rotation curves whose typical pattern after a few 
kpc from the center differs from the Keplerian \j\fr fall-off expected from 
the usual dynamics applied to the electromagnetically-observed matter. 

As a possible solution of this puzzle, the existence of non-baryonic, 
weakly-interacting Cold Dark (in the sense that its existence is indirectly 
inferred only from its gravitational action, not from emitted electromagnetic 
radiation) Matter (CDM) was proposed to reconcile the predictions with the 
observations [3] in the framework of the standard gravitational physics. 

Oppositely, it was postulated that the Newtonian laws of gravitation 
may have to be modified on certain acceleration scales to correctly account 
for the observed anomalous kinematics of such astrophysical systems with- 
out resorting to still undetected exotic forms of matter. One of the most 
phenomenologically successful modification of the inverse-square Newtonian 
law, mainly with respect to spiral galaxies, is the MOdified Newtonian Dy- 
namics (MOND) [H [6] which postulates that for systems experiencing 
total gravitational accelerations A < Aq, with [7] 

A = (1-2 ±0.27) x 1(T 10 m s" 2 , (1) 

the overall gravitational acceleration felt gets modified according to 



VAqGm _ 

A Amond = r. 



More generally, it hold£] 



A=^f^,x^^; (3) 
fj,(x) A 

fi(x) — > 1 for x S> 1, i.e. for large accelerations (with respect to Aq), while 
fj,(x) — > x yielding eq. ([2]) for x <C 1, i.e. for small accelerations. The most 
widely used forms for the interpolating function \i are 

= t~~ ) m (4) 

1 + X 



1 Actually, eq. is exactly valid for isolated mass distributions endowed with partic- 
ular symmetries. 
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Such forms, and also others, as we will see later, can be reduced to the 
following high- acceleration limit (i> 1) 



/j, « 1 — kox 



a 



(6) 



Indeed, eq. dU corresponds to ko = 1, a = 1, while eq. (0) corresponds to 
ko = 1/2 and a = 2. It recently turned out that the simpler form of eq. 
yields better results in fitting the terminal velocity curve of the Milky Way, 
the rotation curve of the standard external galaxy NGC 3198 [8l[ini[II] and 
of a sample of 17 high surface brightness, early-type disc galaxies |12| . Eq. 
([3]) strictly holds for co-planar, spherically and axially symmetric isolated 
mass distributions [13] ; otherwise, for generic mass densities the full modified 
(non-relativistic) Poisson equation^] [9] 



where U is the gravitational potential, G is the Newtonian constant of grav- 
itation and p is the matter density generating U, must be used. 

Attempts to yield a physical foundation to MOND, especially in terms 
of a relativistic covariant theory, can be be found in, e.g., [HI [151 IS]; for 
recent reviews of various aspects of the MOND paradigm, see [XT } fl8 | \19\. 

After setting the theoretical background which we will use in the rest 
of the paper, we will explore the strong MONDian regime in the remote 
(r = 50 — 150 kAU) periphery of the solar system, where the Oort cloud 
[20] . populated by a huge number of small bodies moving along very eccen- 
tric and inclined to the ecliptic orbits, is supposed to exist; for preliminary 
investigations on such a topic, see Ref. [U |2"T] . In particular, we will nu- 
merically investigate the modifications that MOND would induce on the 
Newtonian orbits of a test particle moving in such a region. We will also 
have to consider the subtle External Field Effect (EFE) which, in MOND, 
takes into account the influence of the Galactic field in the dynamics of the 
bodies of the solar system giving rise, in the regions in which x -C 1, to a 
total gravitational acceleration quite different from the scheme consisting of 
the Newtonian one plus some perturbative correction(s) to it. 

For other works on MOND in the inner regions of the solar system, see 



2 By posing A = — Vf7, eq. ([7]) implies that, in general, fj,(A/Ao)A = ANewton + V X h; 
it can be shown [9] that at great distances from an isolated matter distribution with mass 
M it holds fi(A/A )A = -(GM/r 2 )r + 0(r~ 3 ). 




Ref. [a[ia[23l[Ml[25l[26l[271[Ml[32l[33]. 
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2 The External Field Effect in MOND 



In the framework of MOND, the internal dynamics of a gravitating system s 
embedded in a larger one S is affected by the external background field E of 
S even if it is constant and uniform, thus implying a violation of the Strong 
Equivalence Principle: it is the so-called External Field Effect (EFE). In 
the case of the solar system, E would be A cen rj 10~ 10 m s~ 2 because of its 
motion through the Milky Way [4"1 19"1 12T] . 

Clarifying EFE's concrete effects on the orbital dynamics of s is not an 
easy task also because of some possible misunderstandings due to the intrin- 
sic difficulty of the subject. The existence in the inner regions of the solar 
system of its simplest forrrjl, i.e. a constant and uniform vector field E added 
to the standard Newtonian monopole with E ~ 10 -10 m s -2 , has been ruled 
out in |29j by comparing analytical calculation with the estimated correc- 
tions Am to the standard perihelion precessions of some planets determined 
with the EPM ephemerides by E.V. Pitjeva (RAS, IAA) [30]; it turned out 
that the upper bound on an anomalous acceleration with the characteris- 
tics of E is of the order of 10 -14 m s -2 . Such a result has been recently 
confirmed by W.M. Folkner [31], although in a different contex10, with a 
modified version of the latest DE ephemerides by NASA JPL including a 
constant and uniform extra-acceleration with which he fitted long observa- 
tional data records. Milgrom in [32] has put forth a different form for the 
planetary EFE yielding an additional quadrupolar extra-acceleration which 
mimics the action of distant mass; some consequences have been investigated 
in [33] , 

Moving to the remote periphery of the solar system, which is our present 
target, let us define the following quantities 

1 = ^>1, (8) 



r t = J^- = 6.833 kAU, (9) 



.4, 



o 



L = -(f), (10) 



\x \dx 

= Kv), (11) 

L g = Hn). (12) 



' 5 This was, at least, the view about EFE in the solar system of some researchers active 
in the field. 

4 Actually, he considered an extra-acceelration directed radially towards the Sun. 
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In the weak acceleration regime, for 

r > r t rT 1/2 , (13) 

which is fully satisfied in the Oort cloud, the action of EFE is different, so 
that the total barycentric acceleration felt by an Oort comet is^| [9j [32] 

GM ' ^ +V 2 + A' 3,2 ( T f i |. (14) 



^(1 + ^)1/2 Vl + L 



z 



It must be stressed that, contrary to what it could be thought at first sight, 
eq. (iT4|) is not a small correction to the Newtonian monopole ANowton = 
— {GM/r 2 )r to be tackled with the usual perturbative techniques like, e.g., 
the Gauss equations for the variation of the elements. Indeed, eq. (|14ft 
must be treated as a whole because it is the total gravitational pull felt by a 
body in the weak acceleration regime according to MOND; it fully embeds 
EFE which cannot be disentangled. The only consistent alternative to the 
MONDian acceleration o|§ eq. (|14p is the standard Newtonian term. It may 
be interesting to show, from a purely phenomenological point of view, that 
applying eq. (fl"4"j) to the major bodies of the solar system in its planetary 
regions yields absurd results. Indeed, by considering, e.g., the Earth, the 
discrepancy between eq. (|14p and the standard Newtonian monopole turns 
out to be of the order of 2—8 x 10~ 4 m s~ 2 which corresponds to an enormous 
shift of up to 1 x 10 10 m in the heliocentric radial distance of our planet. 

Note that, since the ecliptic longitude and latitude of the Galactic Center 
are about Aqc ~ 180 deg, /3gc ~ — 6 deg [33J, EFE is directed along the X 
axis of the ICRF, i.e. the barycentric frame in which the motion of solar 
system's objects are usually studied. Concerning L, we have 

x 1 

W = TT^ Ll = TTx' (15) 

^=(1+^)1/2 -+ L * = TTx-^ (16) 
^ 2= (1 + ^3/2)2/3 ^ -+ L 3/2 = TT ^. (17) 



5 It comes from eq. J7J by imposing a particular boundary condition, i.e. that for 
r — > oo, V(7 — > —E, and by assuming A/E <C 1 [9]. 

6 Given the scenario considered here, i.e. the solar system embedded in the Milky Way, 
it would be inconsistent to consider the simple form of the MOND acceleration of eq. @. 
Indeed, it does not encompass EFE, holding only for isolated systems or when E/Ao <S 1 
like, e.g., in the case of some satellites of the Milky Way embedded in the external field 
of M31 Andromeda. 
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The value of rj depends on the Galactic field at the solar system's location 



which can be obtained from its centrifugal acceleration 

V 2 

^cen = (18) 

where V is the speed of the Local Standard of Rest (LSR) and R = 8.5 



kpc is the Galactocentric distance. The standard IAU value for the LSR 
speed is V = 220 km s _1 |34| . but recent determinations obtained with the 
Very Long Baseline Array and the Japanese VLBI Exploration of Radio 
Astronomy project yield a higher value: V = 254 ± 16 km s" 1 [35J. Thus, rj 
ranges from 1.5 to 2.3. 

3 Orbits of Oort comets in MOND 
3.1 Ecliptic orbits 

We will, now, consider an Oort comet whose Newtonian orbit covers the 
entire extension of the Oort cloud. It has semimajor axis a = 100 kAU and 
eccentricity e = 0.5, so that its perihelion is 50 kAU and its aphelion is 150 
kAU; for the sake of simplicity, we will assume it lies in the ecliptic plane. 
Its Newtonian orbital period is Pb = 31.6 Myr. We will, first, use /U 3 / 2 . The 
value rj = 2.0, corresponding to V = 254 km s _1 , yields 

fig = 0.82 (19) 
L g = 0.25 (20) 

Figure [Q depicts the numerically integrated Newtonian (dashed blue line) 
and MONDian (dash-dotted red line) orbits for the same initial conditions. 
The integration has been performed backward in time over one (Keplerian) 
orbital period. In Figure [2] we show the case rj = 1.5 (V = 220 km s _1 ) 
yielding 

jig = 0.75 (21) 
L g = 0.34. (22) 

The MOND trajectories are not closed and the points of minimum and 
maximum distances from the Sun do not coincide with the Newtonian ones. 
In general, they do not show a regular pattern. Moreover, the MOND 
paths are much less spatially extended that the Newtonian ones; the overall 
shrinking of the orbit is more marked for the standard value of the LSR 
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Figure 1: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /U3/2, 
77 = 2.0 (V — 254 km s _1 ). The initial conditions are xo = a(l — e),yo — zo = 0,xo = 
0, 2)0 = na« / Zo = 0. The time span of the integration is — Pb < t < 0. 



circular speed (Figure [2|). Finally, the MOND trajectories experience high 
frequency variations during one (Keplerian) orbital period. Such effects 
are particularly notable for highly elliptical Newtonian orbits, as shown by 
Figure [3r Figure H] for e = 0.9 and — Pb < t < 0. In the Newtonian case, 
their perihelia are at 10 kAU, while the aphelia lie at 190 kAU, so that the 
condition of eq. f|13f) is still well fulfilled. The MOND orbits resemble 
involved clews comprised within 60 kAU x 60 kAU. 

At this point, it must be noted that in Figure [T]-Figure [4] we used the 
state vector at the perihelion as initial conditions for the numerical inte- 
grations. We have to check if the behavior described is somewhat related 
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Figure 2: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /U3/2, 
r\ — 1.5 (V — 220 km s _1 ). The initial conditions are xo = a(l — e),yo = zo = 0,xo = 

0, 2)0 = flaw 73 s , io = 0. The time span of the integration is — Pb < t < 0. 



to the peculiar initial conditions chosen. To this aim, in Figure [5]-Figure [6] 
we consider paths starting at the (Newtonian) aphelion for e = 0.5, while 
the case e = 0.9 is depicted in Figure [Tj-Figure Also in this cases we 
integrated the equations of motion backward in time over one (Keplerian) 
orbital period Pb- It can be noted that the MOND trajectories starting at 
the (Newtonian) aphelion are quite different with respect to those starting 
at the (Newtonian) perihelion depicted before. Again, they are not close 
and no regular patterns occur. In this case, the areas swept by the MOND 
paths are larger than in the previous case, and they occupy a large part of 
the Newtonian ones. For e = 0.9 the MOND orbits approximately lie within 
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Figure 3: Numerically integrated orbits of an Oort comet with a — 100 kAU, e = 0.9, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with ^3/2, 
rj — 2.0 (V = 254 km s _1 ). The initial conditions are xo = a(l — e),yo = zq = 0,±o = 

0, yo = naw yz^-, Zo ~ 0. The time span of the integration is — Pb < t < 0. 
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Figure 4: Numerically integrated orbits of an Oort comet with a — 100 kAU, e = 0.9, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with jJ.3/2, 
rj — 1.5 (V — 220 km s" 1 ). The initial conditions are xo = a(l — e),yo = zo = 0,±o = 

0, yo = naJ Zo = 0. The time span of the integration is — Pb < t < 0. 



180 kAU x 60 kAU. 

The case for jj,2 is rather similar, so that we will not depict the related 
figures for saving space. 

Let us, now, examine the case jii. In Figure [9] we show the trajectory 
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Figure 5: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /U3/2, 
77 = 2.0 (V = 254 km s _1 ). The initial conditions are xo = — a(l + e),yo = zo = 0,io = 
0, yo = — na./Yr^, io = 0. The time span of the integration is — Pb < t < 0. 



due to it of the Oort comet with e = 0.5 over Pb for rj = 2.0 which implies 



The case r] = 1.5, yielding 



ix g = 0.67 (23) 
L g = 0.33. (24) 



Hg = 0.60 (25) 
L g = 0.39, (26) 



is shown in Figure [T0l for — Pb < i < 0. The case of highly elliptic orbits (e 
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Figure 6: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /U3/2, 
77 = 1.5 (V = 220 km s _1 ). The initial conditions are xo = — a(l + e),yo = zo = 0,io = 
0, yo = — na./yj^, Zq = 0. The time span of the integration is — Pb < t < 0. 



0.9) is more intricate, as shown by Figure [IT] and Figure [T2j Indeed, the 
MOND paths resemble confuse clews confined within small spatial regions 
shifted with respect to the fi 3 / 2 case. 

Also the trajectories of Figure [9]-Figure H2] start from the (Newtonian) 
perihelia. If, instead, we choose the (Newtonian) aphelia quite different 
paths occur also for this form of /U. Figure [131-Figure 1161 show them. 
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Figure 7: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.9, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /X3/2, 
77 = 2.0 (V = 254 km s" 1 ). The initial conditions are xo = — a(l + e),yo = zo = 0,io = 
0, yo = — naJ jrr, Zo = 0. The time span of the integration is — Pb < t < 0. 



3.2 Nearly polar orbits 

Let us, now consider the case of orbits showing high inclinations / to the 
ecliptic. For practical reasons, here we will only show some cases. In Figure 
[TTl-Figure [191 we depict the sections in the coordinate planes of an orbit 
with a = 66.6 kAU, e = 0.92, / = 81 deg for fi\ and rj = 1.5. The initial 
conditions chosen are arbitrary in the sense that, contrary to Section 13. 1\ 
neither the perihelion nor the aphelion have been used as starting points. 
The case of /i 3 / 2 and r] = 1.5 is illustrated in Figure [20T- Figure 1221 



4 Consequences on the Oort cloud 

The features of the MOND trajectories shown in the previous pictures sug- 
gest that EFE may have consequences on the interaction of the Oort-like 
objects with passing stars |20| by changing their perturbing effects and, 
thus, also altering the number of long-period comets launched into the in- 
ner regions of the solar system and the number of comets left in the cloud 
throughout its history. Indeed, in the standard picture, the comets moving 
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Figure 8: Numerically integrated orbits of an Oort comet with a — 100 kAU, e = 0.9, 
Pb = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /X3/2, 
rj = 1.5 (V = 220 km s _1 ). The initial conditions are xo = — o(l + e),yo = zo = 0,xo = 

0, yo = —nax jrr, io = 0. The time span of the integration is — Pb < t < 0. 



along very (Newtonian) elongated orbits may come relatively close to a star 
of mass M* suffering a change in velocity Av which approximately is |2L)| 

Aw = -i 27 

where is the star's velocity with respect to the Sun and d is the distance of 
closest approach with the Oort object. Moreover, also the perturbing effects 
of the Galactic tides would be altered. In particular, for those particular 
initial configurations yielding highly shrunken paths with respect to the 
Newtonian case the perturbing effects of the passing stars and of the Galactic 
tides may get reduced. Anyway, it is very difficult to realistically predict 
the modifications that the Oort cloud would undergo under the action of 
MOND with EFE. Extensive numerical simulations like, e.g., the one in [36] 
performed in a different context, would be required; they are beyond the 
scope of the present paper. 
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Figure 9: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
P b = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with fj,i, t] = 2.0 
(V = 254 km s _1 ). The initial conditions are xo = a(l — e),yo = zo = 0, xq = 0,yo = 

naJ Y^f , «o = 0. The time span of the integration is — Pb < t < 0. 



5 Summary and conclusions 

We investigated the orbital motions of test particles according to MOND 
(with EFE) in the Oort cloud (r « 50 - 150 kAU). As MONDian interpo- 
lating function fJ,(x), we extensively used the forms \x\ = 1/(1 +x),fi2 = 
x/\/l + x 2 , /Z3/2 = a?/ (1 + x 3//2 ) 2 / 3 . We numerically integrated both the 
MOND and the Newtonian equations of motion in Cartesian coordinates 
sharing the same initial conditions backward in time over one (Keplerian) 
orbital period. We considered both ecliptic and nearly polar trajectories, all 
corresponding to high (Newtonian) eccentricities (e = 0.5 — 0.9). In order to 
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Figure 10: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
Pb = 31.6 Myr. Continuous blue line: Newton. Dash-dotted red line: MOND with /ii, 
rj = 1.5 (V = 220 km s _1 ). The initial conditions are xo = a(l — e),yo = zo = 0,±o = 

0, yo = naJ j^-, Zo = 0. The time span of the integration is — Pb < t < 0. 



evaluate the characteristic MONDian EFE parameters [i g and L g entering 
the problem, we used two different values (V = 220 km s _1 and V = 254 km 
s _1 ) of the circular speed of the solar system's motion through the Galaxy 
yielding the Milky Way's gravitational field at the Sun's location. 

Our results show that the MOND orbits are quite different with respect 
to the Newtonian ones; in general, they are open trajectories which do not 
exhibit any regular patterns. Moreover, they are highly sensitive to the 
initial conditions in the sense that different sets of state vectors yielding 
the same Keplerian ellipses generate quite different MONDian paths. For 
certain initial configurations in the ecliptic, corresponding to Newtonian 
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Figure 11: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.9, 
P b = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with fn, r) = 2.0 
(V — 254 km s" 1 ). The initial conditions are xo = a(l — e),yo = zo = 0,io = 0, yo = 

naJ Y^f , io = 0. The time span of the integration is — Pb < t < 0. 
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Figure 12: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.9, 
P b = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /xi, 77 = 1.5 
(V = 220 km s" 1 ). The initial conditions are xo = a(l — e),yo = zo = 0,io = 0, yo = 

na \/ T^f i = 0. The time span of the integration is — Pb < t < 0. 



perihelion passages of highly eccentric orbits, the MOND trajectories reduce 
to intricate clews spanning very small spatial regions with respect to the 
Newtonian case. For other initial configurations, both in the ecliptic and 
outside it, the MOND trajectories are completely different: among other 
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Figure 13: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
P b = 31.6 Myr. Dashed blue line: Newton. Dash-dotted red line: MOND with /xi, 77 = 2.0 
(V — 254 km s _1 ). The initial conditions are xo = —a(l + e),yo — zo — 0,xo = 0, yo = 

—na,\ jt^, io = 0. The time span of the integration is — Pb < t < 0. 



things, their spatial extensions are larger than in the ecliptic perihelion 
cases. As a consequence, the structure and the dynamical history of the 
Oort cloud, in the deep MONDian regime, would be altered with respect to 
the standard Newtonian picture in a way which is difficult to realistically 
predict. Extensive numerical simulations of the dynamics of the Oort cloud 
in MOND with EFE would be helpful. 
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Figure 14: Numerically integrated orbits of an Oort comet with a = 100 kAU, e = 0.5, 
Pb = 31.6 Myr. Continuous blue line: Newton. Dash-dotted red line: MOND with /xi, 
r\ — 1.5 (V = 220 km s _1 ). The initial conditions are xo = — a(l + e),yo = zo = 0,io = 
0, yo = -naWTv 2 , io = 0. The time span of the integration is — Pb < t < 0. 
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Figure 19: Section in the coordinate {YZ} plane of the numerically integrated orbits of 
an Oort comet with a = 66.6 kAU, e = 0.92, / = 81 deg. Dashed blue line: Newton. Dash- 
dotted red line: MOND with /ii , r\ = 1.5 (V = 220 km s" 1 ). The initial conditions are xo = 
40 kAU,2/ = 30 kAU,z = 5 kAU,i = -23 kAU Myr" 1 ,^ = -15 kAU Myr~\io = 
— 15 kAU Myr -1 . The time span of the integration is —Pb < t < 0. 
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Figure 20: Section in the coordinate {XY} plane of the numerically integrated or- 
bits of an Oort comet with a — 66.6 kAU, e = 0.92, / = 81 deg. Dashed blue line: 
Newton. Dash-dotted red line: MOND with p 3/2 , rj = 1.5 (V = 220 km s" 1 ). The ini- 
tial conditions are x = 40 kAU, yo = 30 kAU, z = 5 kAU,±o = —23 kAU Myr -1 , y = 
-15 kAU Myr -1 , z = -15 kAU Myr -1 . The time span of the integration is ™P b < t < 0. 
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Figure 21: Section in the coordinate {XZ} plane of the numerically integrated orbits 
of an Oort comet with a — 66.6 kAU, e = 0.92, I = 81 deg. Dashed blue line: New- 
ton. Dash-dotted red line: MOND with fj, 3/2 , r] = 1.5 (V = 220 km s -1 ). The initial 
conditions are x = 40 kA\J,y = 30 kAU, z = 5 kAU,±o = —23 kAU Myr _1 ,y = 
-15 kAU Myr -1 , z = -15 kAU Myr -1 . The time span of the integration is -P b < t < 0. 
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Figure 22: Section in the coordinate {YZ} plane of the numerically integrated orbits 
of an Oort comet with a — 66.6 kAU, e = 0.92, I = 81 deg. Dashed blue line: New- 
ton. Dash-dotted red line: MOND with /i 3/2 , T] = 1.5 (V = 220 km s" 1 ). The initial 
conditions are x = 40 kAU,yo = 30 kAU, z = 5 kAU,i = —23 kAU Myr _1 ,y = 
-15 kAU Myr -1 , z = -15 kAU Myr" 1 . The time span of the integration is -P b < t < 0. 
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